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Abstract 

fH We introduce a framework for characterizing Matrix Product States 

O | (MPS) and Projected Entangled Pair States (PEPS) in terms of symme- 

1 1 ^ tries. This allows us to understand how PEPS appear as ground states of 

local Hamiltonians with finitely degenerate ground states and to charac- 
terize the ground state subspace. Subsequently, we apply our framework 
to show how the topological properties of these ground states can be ex- 
O"* plained solely from the symmetry: We prove that ground states are locally 

indistinguishable and can be transformed into each other by acting on a 
restricted region, we explain the origin of the topological entropy, and we 
discuss how to renormalize these states based on their symmetries. Fi- 
nally, we show how the anyonic character of excitations can be understood 
as a consequence of the underlying symmetries. 

oo 
m 

^h" 1 Introduction 

O 

What are the entanglement properties of quantum many-body states which 
characterize ground states of Hamiltonians with local interactions? The answer 
seems to be "an area law" : the bipartite entanglement between any region and 
its complement grows as the area separating them - and not as their volume, as 
is the case for a random state (see [l] for a recent review). Moreover, particular 
corrections to this scaling law are linked with critical points (logarithmic cor- 
rections) or topological order (additive corrections). A rigorous general proof 
of the area law, however, could up to now only be given for the case of one- 
dimensional systems [2] , where an area law has been proven for all systems with 
an energy gap above the ground state, whereas the currently strongest result 
for two dimensions [3j|4j requires a hypothesis on the eigenvalue distribution 
of the Hamiltonian. Surprisingly, there is a completely general proof in arbi- 
trary dimensions if instead, we consider the corresponding quantity for thermal 
states [5], and similar links to topological order persist [6j. 
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The area law can be taken as a guideline for designing classes of quantum 
states which allow to faithfully approximate ground states of local Hamiltoni- 
ans. There are several of these classes in the literature: Matrix Product States 
(MPS) |7 and Projected Entangled Pair States (PEPS) [I] are most directly 
motivated by the area law, but there are other approaches such as MERA (the 
Multi-Scale Entanglement Renormalization Ansatz) [9] which e.g. is based on 
the scale invariance of critical systems; all these classes are summarized under 
the name of Tensor Network or Tensor Product States. Though the main mo- 
tivation to introduce them was numerical - they consitute variational ansatzes 
over which one minimizes the energy of a target Hamiltonians and thus obtains 
an approximate description of the ground state - they have turned out to be 
powerful tools for characterizing the role of entanglement in quantum many 
body systems, and thus helped to improve our understanding of their physics. 

In this paper, we are going to present a theoretical framework which al- 
lows us to understand how MPS and PEPS appear as ground states of local 
Hamiltonians, and to characterize the properties of their ground state subspace. 
This encompasses previously known results for MPS and particular instances 
of PEPS, while simultaneously giving rise to a range of new phenomena, in 
particular topological effects. Our work is motivated by the contrast between 
the rather complete understanding in one and the rather sparse picture in two 
dimensions, and we will review what is known in the following. We will thereby 
focus on analytical results, and refer the reader interested in numerical aspects 
to 
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1.1 Matrix Product States 

Matrix Product States (MPS) [7] form a family of one-dimensional quantum 
states whose description is inherently local, in the sense that the degree to 
which two spins can be correlated is related to their distance. The total amount 
of correlations across any cut is controlled by a parameter called the bond dimen- 
sion, such that increasing the bond dimension allows to grow the set of states 
described. MPS have a long history, which was renewed in 1992 when two appar- 



ently independent papers appeared: In 11 , Fannes, Nachtergaele, and Werner 



generalized the AKLT construction of 12 by introducing the so-called Finitely 



Correlated States, which in retrospect can be interpreted as MPS defined on 
an infinite chain; in fact, this work layed the basis for our understanding of 
MPS and introduced many techniques which later proved useful in characteriz- 
ing MPS [13]. The other was [l4j, where White introduced the Density Matrix 
Renormalization Group (DMRG) algorithm, which can now be understood as 
a variational algorithm over the set of MPS. In [15], MPS were explained from 
a quantum information point of view by distributing "virtual" maximally en- 
tangled pairs between adjacent sites which can only be partially accessed by 
acting on the physical system. This entanglement-based perspective has since 
then fostered a wide variety of results. 
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I. The complexity of simulating ID systems 

Motivated by the extreme success of DMRG, people investigated how hard or 
easy the problem of approximating the ground state of a ID local Hamiltonian 
(or simply its energy) was. The history of this problem is full of interesting 
positive and negative results. A number of them was devoted to prove that 
every ground state of a gapped ID local Hamiltonian can be approximated 



by an MPS 16 17 ; this was finally proven by Hastings [2j, justifying the use 
of MPS as the appropriate representation of the state of one-dimensional spin 
systems. Very recently, also in the positive, it was shown that dynamical pro- 
gramming could be used to find the best approximation to the ground state of 
a one-dimensional system within the set of MPS with fixed bond dimension in 
a provably efficient way 18 19 . On the other hand, in the negative it could be 



shown that finding the ground state energy of Hamiltonians whose ground states 
are MPS with a bond dimension polynomial in the system size is NP-hard 20 ; 
this is based on a previous result of Aharonov et al. (2ll proving that finding 
the ground state energy of ID Hamiltonians is QMA-complete (the quantum 
version of NP-complete) . 



II. Hidden orders, symmetries and entanglement in spin chains 

As we have seen, MPS provide the right description for one-dimensional quan- 
tum spin chains. Therefore, and given their simple structure, one can employ 
MPS to improve our understanding of the physics of one-dimensional systems. 
One field in which significant insight could be gained was the characterization 
of symmetries in terms of entanglement. First, the relation between string or- 
der parameters and localizable entanglement was explained in 22 23 . In 24 



(see also [25]), global symmetries in generic MPS have been characterized, and 
related to the existence of string order parameters, thus explaining many of the 
properties of string order, for instance its fragility 
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of global symmetries was generalized to arbitrary MPS in 27 



This characterization 
where it was 



used to shed light on the Hamiltonian-frce nature of the Lieb-Schultz-Mattis 
theorem as well as to find new SU(2)-invariant two-body Hamiltonians with 
MPS ground states, beyond the AKLT and Majumdar-Ghosh models. Other 



examples of MPS with global symmetries were already provided in 11.28 29 



Recently, also reflection symmetry has been investigated, showing how it pro- 
vides topological protection of some MPS such as the odd-spin AKLT model, 
as opposed to the even-spin case |30] . 

MPS have also been extremely useful in understanding the scaling of entan- 
glement in quantum spin chains, where special attention has been devoted to 



the case of quantum phase transitions. In 31 , MPS were used to give examples 



of phase transitions with unexpected properties, namely analytic ground state 
energy and finite entanglement entropy of an infinite half-chain; the entangle- 
ment properties of these examples were further analyzed in 32033 .In 34 , MPS 



theory was used to compute how the geometric entanglement with respect to 
large blocks diverges logarithmically with the correlation length near a critical 
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point, and thus takes the role of an order parameter (see also [33] ) . 

Apart from that, MPS theory has been used to decompose global operations 
(such as cloning or the creation of an entangled state) into a sequence of local 



operations 35 36 , to characterize renormalization group transformations and 



their fixed points in ID 37 , to understand which quantum circuits can be 



simulated classically 38 , or even to propose new numerical methods to solve 



differential equations 39 or to compress images 40 



But what is it that makes MPS so useful in deriving all these results? 
III. The structure of MPS 

The main reason seems to be that MPS, despite being able to faithfully represent 
the states of one-dimensional systems, have a simple and well-understood struc- 
ture, which makes them quite easy to deal with. For instance, as shown in |38] , 
they naturally reflect the Schmidt decomposition at any cut across the chain, 
which makes dealing with their entanglement properties particularly easy. If 
one moreover restricts to the physically relevant case of translational invariant 
states, it turns out that one can fully characterize the set of all translation- 



ally invariant MPS by bringing them into a canonical form 11 13 ; in fact this 
canonical form constitutes one of the main ingredients in many of the results 
mentioned above. 

In the canonical form, the matrices characterizing the MPS obtain a block 
diagonal form, and the properties of the state can be simply read off the struc- 
ture of these blocks. In particular, for the case of one block (termed the in- 
jective case), it can be shown that the MPS arises as the unique ground state 
of a so-called "parent" Hamiltonian with local interactions, which moreover is 
frustration free. For the non-injective case where one has several blocks, the 
number of blocks determines the degeneracy of the parent Hamiltonian, and 
the ground state subspace is spanned by the injective MPS described by the 



individual blocks 13 . Beyond that, the injective case has other nice properties, 



such as an exponential decay of correlations. For the case of an infinite chain, 



all these properties were proven in 11 , together with the fact that the parent 



Hamiltonian of an injective MPS has an energy gap above the ground state. 
The block structure of the canonical form is also useful beyond the relation of 
Hamiltonians and ground states, and e.g. allows to read off the type of the RG 
fixed point of a given ID system. 

1.2 PEPS 

Projected Entangled Pair States (PEPS) constitute the natural generalization 
of MPS to two and higher dimensions, motivated by the quantum information 
perspective on MPS which views them as arising from virtual entangled pairs 
between nearest neighbors 1 8 , 41 . Though there has not yet been a complete 



formal proof that PEPS approximate efficiently all ground states of gapped local 
Hamiltonians, this could be proven under a (realistic) assumption on the spectral 



density in the low-energy regime 42 , showing that PEPS are the appropriate 
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class to describe a large variety of two-dimensional systems. However, as com- 
pared to MPS, PEPS are much harder to deal with: For instance, computing- 
expectation values of local observables, which would be the key ingredient in 
any variational algorithm such as DMRG, is a #P-complete problem, and thus 
in particular NP-hard [43]. This poses an obstacle to numerical methods, and 
different ideas to overcome this problem have been proposed (we refer again 
to 10 for numerical issues). Fortunately, the bad news comes with good ones: 
The increase in complexity allows to find a much larger variety of different in- 
teresting behavior within PEPS as compared to MPS; for instance, in [44] it is 
shown that there exist PEPS with a power-law decay of two-point correlation 
functions, something which cannot be achieved for MPS. 

I. Many examples 

There have been identified various different classes of PEPS which exhibit rich 
properties. To start with, in [45] it has been shown how the interpretation 
of the 2D cluster state as a PEPS can be used to understand measurement 



based quantum computation 46 - a way of performing a quantum computation 
solely by measurements - by viewing it as a way to carry out the computation 
as teleportation-based computation on the virtal maximally entangled states 
underlying the PEPS description. This motivated the search for different models 
for measurement based quantum computation, and indeed models with very 



different properties have been subsequently proposed 47-49 . 

Another category of examples has been found with regard to topological 
models, where it was realized that many topologically ordered states have a 
PEPS representation with a small bond dimension: The case of Kitaev's toric 



code was observed in 44 , and this was later generalized to all string net models 
in 50 . Yet, despite their ability to describe these states, up to now PEPS did 
not help much in understanding the topological behavior of these states, which 
is one of the things we will assess in this work. 

II. Few general results 

As we have seen, the class of PEPS is rich enough to incorporate states with 
a variety of different behaviors. However, given the complexity of topological 
systems or of (measurement based) quantum computation, is this class still 
simple enough to prove useful as a tool to improve our understanding of 2D 
systems? That is, can PEPS help to uncover new effects and relations in nature, 
or to give a better understanding of the mechanisms behind quantum effects in 
two dimensions? 

Judging from the experience with one-dimensional systems, in order to do so 
it would be highly desirable to have an understanding of the structure of PEPS 
comparable to the one obtained using the canonical form in one dimension. 
As it turns out, in the case of "injective" PEPS, several ID results can be 
transferred; in particular, injective PEPS appear as unique ground states of their 



parent Hamiltonian 51 . Also, in 52 it is shown how global symmetries can 
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be characterized in injective PEPS, which helped to understand the mechanism 



behind the 2D version of the Lieb-Schultz-Mattis theorem 53 , to define an 



appropriate analogue for string orders in 2D 24 , and to improve the PEPS- 
based algorithms used to simulate 2D systems with symmetries 54 . While 
these results illustrate that PEPS are a useful tool to understand properties of 
quantum states which appear as unique ground states of local Hamiltonians, it is 
also true that some of the most interesting physics in two dimensions takes place 
in systems which do not have unique ground states, but rather lowly degenerate 
ones, such as systems with symmetry broken phases or states with topological 
order. 

Yet, in order to be able to fully apply the toolbox of PEPS to the understand- 
ing of these systems, it would be crucial to have a mathematical characterization 
of the structure of PEPS, providing a framework similar to the one which proved 
so useful for one-dimensional systems: Is there a canonical form for PEPS which 
allows to easily determine their properties, and how can it be found? How do 
PEPS appear as ground states of local "parent" Hamiltonians, and what is the 
ground state degeneracy? What is the structure of the ground state subspace, 
and how do these states relate to the PEPS under consideration? 



1.3 Content of the paper 

In this paper, we introduce a framework for characterizing both MPS and PEPS 
in terms of symmetries of the underlying tensors. This classification allows us to 
rederive all results known for one-dimensional systems, while it can be equally 
applied to the characterization of two- and higher-dimensional PEPS states, 
answering the aforementioned questions: Using the characterization based on 
symmetries, we prove how PEPS appear as ground states of local Hamiltonians 
with finitely degenerate ground states, and how these states can be obtained as 
variants of the original state. Subsequently, we demonstrate the power of our 
framework by using it to explain in a simple and coherent way the topological 
properties of these ground states: We prove that these states are locally indis- 
tinguishable and can be transformed into each other by acting on a restricted 
region, we explain the origin of the topological entropy, and we discuss how to 
rcnormalize these states based on their symmetries. We also discuss the excita- 
tions of these Hamiltonians, and demonstrate how to understand their anyonic 
statistics as a consequence of their symmetries. Thus, the characterization of 
PEPS in terms of symmetries provides a powerful framework which allows to 
explain a large range of their properties in a coherent and natural way. 

The material in this paper is structured as follows: In Sec. [2j we introduce 
MPS and PEPS and discuss their basic properties. In Sec. [3j we review the 
proof for the "injective" case, in which the MPS appears as a unique ground 
state of the parent Hamiltonian. In Sec. [4j we show how general MPS can be 
classified in terms of symmetries, and use this description to generalize the re- 
sults of the preceding section; we also discuss how our results relate to the ones 
obtained using the original canonical form of [13) . In Sec. [5] we generalize the 
symmetry-based classification of MPS to the case of PEPS, where we derive the 
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parent Hamiltonian and characterize the structure of its ground state space. In 
Sec. [6j we consider a more restricted class of PEPS with symmetries, for which 
we derive a variety of results concerning the structure of the ground state space 
and the parent Hamiltonian, such as topological entropy, local indistinguisha- 
bility, renormalization transformations, or the fact that the parent Hamiltonian 
commutes. We conclude the section by characterizing the excitations of the 
system and explaining how their anyonic statistics emerges from the symme- 
tries of the PEPS. Finally, in Sec. [7J we discuss examples which illustrates the 
applicability of our classification, before we conclude in Sec.[8j 



2 MPS and PEPS 

We start by defining Matrix Product States, which describe the state of a one- 
dimensional chain of d-level systems of length L. 

Definition 2.1. A state \ip) € (C d )® L is called a Matrix Product State (MPS) 
if it can be written as 



\M(A))= J2 tr^.-.^lix, 

lit 

Here, A 1 6 L(C D ) (the space of D x D matrices ot 
bond dimension. 



,*L> • (2-1) 
C ), and D is called the 



Note that we restrict to the translational invariant setting with only one 
tensor A, whereas a general MPS can be defined with a different tensor at 
each site. However, for each MPS describing a translationally invariant state 



a translationally invariant description (2.1| can be constructed; note, however, 
that the bond dimension can increase with L. 

In a more graphical representation, A can be written as a three-index tensor, 



a=J2(a% 



ia/3 



Here, p refers to the "physical" index [characterizing the actual state in (2.1 1] 



and v to the "virtual" index, which is only used to construct the MPS (2.1 1 and 
does not appear in the final state. For clarity, we assign arrows to the virtual 
indices, pointing from bras (-| to kets |-). Connecting the legs of two tensors 
(called a "bond") denotes contraction, 

f o , 

5>%„(^), 7 \i,j) p \a) v (7L = WW 



Thus, we can write the MPS (|2.1|) as 

Ha) = £ 



fl 


h 


h 


h 


A 




A 




A 




A 
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Having this picture in mind, we can immediately define a two-dimensional 
generalization: 

Definition 2.2. A state \ip) € (C d )®( ixL ) is called a Projected Entangled Pairs 
State (PEPS) of bond dimension D if it can be written as 



(2.2) 



with A % a p lS , i — 1, . . . , d, a, /3,7, 5 = 1, . . . , D a five-index tensor, where the 
index denoted i corresponds to the physical system and "goes out of the paper". 



Note that in the definition (2.2 1 we have implicitly defined the direction of 



the arrows, and thus the assignments of kets and bras in A: 



A = J2Ai trb \i) p \l,t) v {r,b\ v = l^i; 



iltrb 



lb 



Let us now introduce two simplifications to the MPS tensors, which will 
allow us to characterize all MPS based on their non-local properties. 

First, we define the "projector" - the map which creates the physical system 
from the virtual layer. 



Definition 2.3. For an MPS given by a tensor A, 

Ml 



V{A) :=$>i 

ict(3 



/SI*V 



(2.3) 



is the map which maps the virtual system to the physical one. The definition 
extends directly to PEPS. 

The importance of the map V(A) lies in the fact that it tells us how virtual 
and the physical space of the MPS are connected. On the one hand, it tells us 
which physical states can be created by acting on the virtual degrees of freedom, 
and on the other hand, it tells us which virtual configurations - which in turn 
enforce physical configurations on the surrounding sites - can be realized by 
acting on the physical system. Thus, studying the properties of V(A) will allow 
us to infer properties of the MPS. 

Let us now see which simplifications we can make in the analysis of V(A), 
given that we are only interested in the non-local properties of the MPS (PEPS). 
Define 71(A) ;= rgV(A), and V(A) := (kerV(A)) 1 - = span({A l } 4 ). Then, 7>(A) 
can be inverted on T>(A) and 71(A), respectively: 



BTiA)- 1 : V(A)~ 1 P(A) 



1\ V{A) A V{A)V(A)- 1 



\K(A) 
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Note that TZ(A) characterizes the local support of the MPS: Projecting sites 
2, . . . , L on any basis state ■ ■ ■ , %l leaves us with J2i 1 ti[A ll X] \ii) (where 
X = A 12 ■ ■ ■ A lL ), and thus, the single-site reduced operator is supported on 
1Z(A). Thus, we can restrict the MPS to the subspace 7Z(A)® L , i.e., each local 
system can be mapped to a dim7?.(^4)-dimensional system by a local isometry, 
i.e., without changing any non-local properties. This implies that we can w.l.o.g. 
restrict our analysis to MPS with the following property. 

Observation 2.4. Any MPS/PEPS \M{A)) can be characterized (up to local 
isometries) by a tensor A for which 'P(A) has a right inverse (denoted by A~ x 
in the diagram) such that: 



A 


]- 


and 


A 1 


A 1 






A 



(2.4) 



where Ht>(A) = ^\x>(A) * s ^ e orthogonal projector on D{A). W.l.o.g., we will 
assume a description with this property from now on. 

The purpose of this paper is to characterize MPS and PEPS by looking at 
the structure of the subspace T>(A), and especially at its symmetries. As it 
turns out, it is sufficient to consider two cases: First, T>(A) — span({A l }i) = 
L(C D ), the space of all linear operators on C D , and second, the case where 
T>(A) = span({A l }j) is an arbitrary C*-algebra, i.e., a linear space of matrices 
closed under multiplication and hermitian conjugation. 

In order to see why we only need to consider these two cases, take an MPS 
with tensor A, and group its sites into super-blocks of k sites each. This results 
in a new MPS with tensor 



B = A A 



(2.5) 

The map V{B) for this new tensor goes from a D 2 to a d fe -dimensional space: 
Thus, for k > 2 log d D, we have that D 2 < d k , which means that the map V(B) 
will typically be injective, and thus D(A) — L(C D )^ These MPS are known 
as injective, and they appear as unique ground states of their associated parent 
Hamiltonians. 

The second case - D (A) being a C*-algebra - arises e.g. if the A 1 are block 
diagonal matrices, but within each block span the whole space of linear opera- 
tors. As shown in [13], this does in fact cover the case of a general MPS, as any 
MPS can be brought into this form. While these states are no longer unique 
ground states of local Hamiltonians, their parent Hamiltonians have a finite 
ground state degeneracy, and the ground states are described by the individual 
blocks of the Ai. 

In the following, we will first review the situation of injective MPS, and 
show how to prove that they are unique ground states of local Hamiltonians. 



1 See [13] for a discussion of how to understand "typical" in this context. 
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We will then turn towards the case where T>(A) is a general C*-algebra, which 
we translate into a condition on the symmetry of A under unitaries. While in 
one dimension this reprodu ces the results previously derived using the block 
structure of the matrices A % 
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13 



it will enable us to generalize these results 
to the two-dimensional scenario, where we will find states exhibiting topological 
order and anyonic excitations, all of which can be understand purely in terms 
of symmetries. 



3 MPS: the injective case 
3.1 Definition and basic properties 

We start by analyzing the injective case in which D(A) — span({A l }i) = L(C D ). 
According to Observation 2.4 we can choose A such that V(A) is invertible. 



This leads us to the following formal definition of injective. 

Definition 3.1 (Injectivity). A tensor A is called injective ifV(A) has a left 



rq -- 
" A - 



(3.1) 



(The corresponding MPS \A4(A)) will also be termed injective.) Intuitively, 
injectivity means that we can achieve any action on the virtual indices by acting 
on the physical spins. 



Lemma 3.2 (Stability under concatenation). The injectivity property (3.1) is 
stable under concatenation of tensors: If A and B are injective, then the tensor 
obtained by concatenating A and B is also injective, since 



A B 



B 



Note that we generally omit normalization constants in the diagrams (the 
contraction of the loop is tr 1=1?). 



3.2 Parent Hamiltonians 

Let us now see how injective MPS give rise to parent Hamiltonians, to which 
they are unique ground states. To this end, note that the two-particle reduced 
operator of the MPS \M{A)) is given by 



P [2] (A) = 



with ip 



[2] 



(A) the projector onto the two-particle state obtained by projecting 



sites 3, . . . , L to the basis state 13, . 



,1L, 



(A)) = J2tT[A^A i2 X is '-^]\i u i 2 ) 



10 



with X 1 



= A 13 ■ ■ ■ A lL G D{A). Thus, p^(A) is supported on the subspace 



5 3 = ^tr[A 4 ^X]|t > j; 



lelf 




(3.2) 



Moreover, the injectivity of A implies that X %3 ''">' lL spans the space of all D x D 
matrices, and thus, p' 2 ' has actually full rank on <S 2 - Analogously to (3.2), one 
can define a sequence of subspaces 



S k = { ]T tx[A i ^--A i "X}\i 1 ,. 



X G L 



(C D )} 



which by the same arguments exactly support the fc-body reduced operator 

The idea for obtaining a parent Hamiltonian is now as follows: Define a two- 
body Hamiltonian which has 52 as its ground state subspace, and let the parent 
Hamiltonian be the sum of these local terms. The proof consists of two parts: 
First, we show that for a chain of length fc with open boundaries, the ground 
state subspace is Sk (i.e., optimal), and second, when closing the boundaries, 
the only state remaining is \A4(A)). This is formalized in the following two 
theorems. 

Theorem 3.3 (Intersection property) . Let A and B be injective tensors. Then, 



A B 



M 



M \ n 



N 



B A 




N > = 



X 



with M,N,X G L(C D ) 



(3.3) 



Proof. It is clear that the right side is contained in the intersection, since for 
any X we can choose 



M 



[A] 

fx} 



and 



N = 



(3.4) 



Conversely, for any \ip) in the l.h.s. of (3.3 1, there exist M, N such that 



\i>) 



A B 



M 



N 



1 L X 

UB A 



Applying the left inverse of V(A) and V(B) to the i and j index, we find that 



N 



N 







B 


-A 



bW^ 



A B 



M 



M 



i.e., N is of the form ( 3.4 1, and thus is contained in the r.h.s. of (3.3). □ 
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Theorem 3.4 (Closure property). For injective A and B, 



B 



M 



M \ n 




N\ = 



B 



(3.5) 



Proof. As in the proof of Theorem |3.3[ the r.h.s. is trivially contained in the 
intersection by choosing N — M — 1. Conversely, by taking an arbitrary 
element in the intersection and applying the left inverses of V(A) and V(B), we 
find that 



M 







A* 


J L 


5 H 




A* 


B 




A 


] 




B 




A 


M 






N 



N 



which proves (3.5 1 



□ 



Let us now put Theorems 3.3 and 3.4 together to show that the MPS |jV1(A)) 



arises as the unique frustration free ground state of a local Hamiltonian. 



Theorem 3.5 (Parent Hamiltonians) . Let A be injective, andS2 as in Eq. (3.2). 

Define 

hi = if (i_1) ® (i - n s J ® \\ 



I L-i-2 



as the orthogonal projector on the subspace orthogonal to S 2 on sites i and i + 1 
(modulo L). Then, 

L 

-^par — ^ ^ hi 



has \M(A)) as its unique and frustration free ground state. 

Proof. For k < L, define j[k - 1] := (i 2) . . . , i k _ x ) and B^ k ~^ := A 12 ■ ■ ■ A 1 "- 
Rewriting 



S, 



k-l 



{ J2 ^[A^B^-^M^Wh^lk-llh 

ilj[k-l],ik 



M G 



(and similarly for C d <E) Sf.-i and Sk), Theorem 3.3 implies that Sk-i <8> C d H 
C d ® Sk-i = Sk, and thus by induction 



Sl — So 



: d )®( fc - 2 ) n c d ® s 2 ® (c d )®( fc - 3 ) n •■■ n (c d )® (fe ~ 2) ®s 2 , (3.6) 



i.e., the subspace supporting the length L chain is given by the intersection of 
the two-body supports S 2 - 
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Now let iTieft = hi + \h 2 A \-\h L -i, and Fright = \h 2 H h \h L -i + h L 

As the hi are projectors, the null space of i?ieft is given by the intersection (3.6) 
i.e., by 



Sieft =Sl = { J2 ^ h B j[k] M] \h,j[k}) 



Correspondingly, the null space of -ff r ; g ht is 

5 right = { tT[A^NB^]\ii,j[ 
hJ[k] 



N e L 



(C D )} , 



and thus by the closure property, Theorem 3.4 |.M(A)) = 6>i c f t n 6> r i g ht is the 
unique zero-energy (i.e., frustration free) ground state of -ff par = H\ c f t + -Bright- 

□ 



Note that all these results hold equally for injective PEPS 51 . We will 



discuss the case of PEPS in detail for the non-injective scenario, which includes 
the injective one as a special case. Note also that the results of this section do 
not rely on the translational invariance of the PEPS - the central Theorems [3^3 
and|3.4| hold for any pair ^4, B of injective tensors. 



4 MPS: the G-injective case 
4.1 Definition and basic properties 

In the following, we will consider the case where A is not injective, but where 
nevertheless 

V = V{A) = [Y J ^ Al \ X ^ £ ] 

i 

forms a C* -algebra. 

In the following, we will characterize the structure of T>, and thus of A, in 
terms of symmetries. Using Observation 2.4 - that V(A) has a left inverse on 
T> - together with an appropriate characterization of T> will allow us to base 
proofs on this left inverse, similar to the injective case. 

Theorem 4.1. For any C* -algebra T> C L(C D ) there exists a finite group G 
and a unitary representation g n- U g such that T> is the commutant ofU g , i.e., 

V = {X e L(C D )\[X 7 U g ] = V 5 € G} . (4.2) 

Proof. Any C* -algebra T> can be decomposed as 

D^l^LfC" 1 '), (4.3) 
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where = denotes unitary equivalence, A = B 3 W, WW^ = 1 : A = WBW^ . 
Now choose a finite group G and a representation U g with irreducible represen- 
tations D l of dimensions di and multiplicity m%, respectively^] 



Now Schur's lemma implies Eq. (4.2 1. 



□ 



Note that conversely, any unitary representation U g has a C*-algebra as 



its commutant (4.2), since f/j 



Thus, the characterization in terms of 



unitaries is equivalent to the characterization in terms of T>, the span of the A 1 . 
This motivates the following definition. 

Definition 4.2. Let g i— > U g be a unitary representation of a finite group G. 
We say that an MPS tensor A is G-injective if 



i) Vi,g: U g A i U^ = A i , 

ii) the map V(A) [EqM 



of U g -invariant matrices, 



has a left inverse on the subspace 
S = {X\[X,U g }=0} 

V-\A)V{A) = l\ s . 



(4.4) 



(4.5) 



In graphical notation, we will denote unitary representations as circles la- 
belled U g , or simply g when unambigous. Note that the arrow now point towards 
the ket on which to apply U g . Condition i) of Definition 



4.2 



then reads 



Note that the group G is more important than its representation U g (though 
we will require certain properties at some point), and in fact, a different repre- 
sentation can be attached to each bond. 



Lemma 4.3. The orthogonal projector cr onto the U g -invariant subspace (4.4) 
is given by 



Here, \G\ is the cardinality of the group G. Thus, condition (4.5) corresponds 
to 

HZ 



-A 



= L 

8 



(We generally omit normalization in diagrams.) 



2 This can be achieved, e.g., by choosing / finite groups G,: which have a di-dimensional 
irreducible representation Di, and letting G = Gi X • • • X Gj, and D'(gi, . . . ,gj) = Di(g{). 
In particular, the choice of G and U g is not unique. 
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Proof. Since 



1 1 gh 

= i^ia E = iif E = < x ) . 

1 1 /ig 1 1 k 

cr is a projection. As it leaves 5 invariant, 

U h tr(X)ul = ^ E ^< = CT W 

(i.e., the image is contained in 5), and it is hermitian, it is the orthogonal 
projector on S. □ 



We will now show the analogue of Lemma 3.2 G-injectivity is stable under 
concatenation of tensors. To this end, we will use the following identity. 

Lemma 4.4. For any unitary representation U g of a finite group, 
with 

1 1 i 

where (di, mi) are dimensions and multiplicites of the irreducible representations 
D l of U g , and A is diagonal in the basis in which U g = @ t D l (g) ® l mi . In a 
formula, 

ti[UlAU g ]U h = J2 tr[U gh -iA]U h = U g . (4.7) 

h h 

Proof. The group orthogonality theorem implies 

g g,m,k,l 

= E 5 i^m,kS m ,i \k) (l\ (4.8) 

1 m,k,l 



\g\ 

di 
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Thus, 

^2tT[U gh -iA]U h =^tT[U k -,A]U k U g 



E 



^TO;tr[L> 4 (£f 



i\G\- 



[Di(k)Di(g)] 



0^(5) 



3 • 



□ 



For a restricted class of representations, we can make a stronger statement. 

Definition 4.5 (Semi-regular representations). A unitary representation g t— > 
U g of a finite group G is called semi-regular if U g contains all irreducible rep- 
resentations ofG. 

Lemma 4.6 (Linear independence of semi-regular representations). Let g h-> U g 
be a semi-regular representation of a group G. Then, 

trplUhA] = 5 gih 



with A as of (4.6). Note that for the regular representation, A oc 1. 

Proof. This follows as |G|tr[[/ fe A] = |G| ^ m^^^D 1 (k)] is the character of 
the regular representation, which is \G\ 5k,i- □ 



We are now ready to prove the analogue of Lemma 3.2 for the G-injective 



case. 



Lemma 4.7 (Stability under concatenation). Let A and B be G-injective ten- 
sors. Then C lJ = A l B J is also G-injective with left inverse 



C ~:=~A" 



B 



(4.9) 



Here, A is defined as in Lemma \4.4\ 
Proof. G-invariance of G follows from 



hr A 



A <hy =)<sy A- A A (sKs) A A - A ~ 



(4.10) 



Moreover, (4.9) is the left inverse of V(C) on the /7 g -invariant subspace since 



A H 




B A 


A 




B 




from Lemma ]AA 



□ 
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4.2 Parent Hamiltonians 



Let us now proceed to the relation of G-injective MPS and parent Hamiltonians. 
The construction is exactly analogous to the injcctivc case: We define 



S k = { E tr[A^.--A i "X}\i 1 



,ik) 



X G U 



and prove the analogues of Theorem 3.3 (Intersection Property) and Theo- 
rem 3.4 (Closure Property), but now for G-injective MPS. While the intersec- 



tion property will be the same as in the injective case, the closure will give rise 
to a subspace whose dimension equals the number of conjugacy classes of G, for 
a properly chosen group G. 



Theorem 4.8 (Intersection property) . Let A, B be G-injective. Then, 



A ~ B 



M 



m \ n 



N 



B~A 




X 



(4.11) 



Proof. Using the G-invariance of A and B, we can infer that we can restrict M, 
N, and X to also be G invariant in the virtual indices. For instance, for any M, 



A'B~ 



M 



this is, we can replace M by 




M , 



M 



the same holds true for N and X. We will always assume symmetrized tensors 
from now on. 

We are now ready to prove (4.11 1. First, the r.h.s. is clearly contained in the 
l.h.s., by choosing AP — A l X, N l = XA 1 . On the other hand, each element in 
the intersection can be simultaneously characterized by a pair N, M of tensors, 
and by applying the inverse maps we find that 



N 



N 















N 




Si 






M 









4>Wl 



M 8 



X M 




M 



This shows that M (and equally N) are of the form required by Eq. (4.11 ), and 
thus proves the theorem. □ 
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Theorem 4.9 (Closure property). For G-injective A and B, 



B 




A 


M 




B 


If 


A 


N 



N 



B~A 



(4.12) 



Before proving the theorem, let us give an intuition why the closure can be 
done using any U g (and nothing else). To this end, regard B as the (L — l)-fold 
blocking of A's. The closures U g are exactly the operators which commute with 
A, and therefore, they can be moved to any position in the chain. Thus, no 
local block of A's needs to hold the closing U g , i.e., the state looks the same 
locally independent of the closure. 

Proof. We may again assume that M and N are G-invariant. It is clear that 
the r.h.s. is contained in the l.h.s., by moving U g to the relevant link and setting 
M = U g or N = U g , respectively. On the other hand, any element in the 
intersection can be written using some M and N, for which it holds that 



M 



zc 
A 



M 



zz 
A 



N 




(In the first step, we have used that M is G-invariant.) Thus, any element of 
the intersection is of the form of the r.h.s. of (4.12), with X g = tr[U g -iNA], 
g = khr x . □ 

Let us now formally define translational invariant MPS with an operator in 
the closure, as they appear in the above theorem. 

Definition 4.10. For an MPS tensor A and K £ L(C D ), we define 

\M(A\K)):= ^ ■■■A^K]\i 1 ,...,i L ) (4.13) 



to be the MPS given by A with closure K 



Theorem 4.11 (Parent Hamiltonians). Let A be G-injective, S2 as in (3.2 1, 

and let 



(1-Hs 2 )®l£ 



The 



L-i-2 



3 Note that one can use a different closure to reduce dramatically the bond dimension. One 
example is the W-state, whose bond dimension grows to infinity with L with the standard 
closure [13||55| , but has bond dimension 2 with a different closure. 
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has a subspace of frustration free ground states spanned by the MPS \M.(A\U g )} 



with U g -closed boundaries. 



Proof. The proof is exactly the same as for Theorem 3.5 except that is is now 
based on Theorems |4.8| and |4.9| the latter leading to the degeneracy of the 
ground state subspace. □ 

Theorem 4.12 (Structure of ground state subspace). Let D 1 ,...,!) 7 be the 

irreducible representations of U g , of dimension <ij. Then, the ground state 
subspace of Theorem is I -fold degenerate, and it is spanned by the MPS 
\M(A\U Z )), where 

^i = ^J2 tT i Di (9' 1 )P g (4.14) 
' ' g 

is the projector onto the subspace supporting the irreducible representation D l 
in U g [proven in (4.8)/. 

Moreover, if U g is a semi-regular representation, I is equal to the number 
of conjugacy classes of G, and the subspace is spanned by the linearly indepen- 
dent states \ A4(A\U g )) , where for each conjugacy class g, one representative is 
chosen. 

Note that the first part of the theorem corresponds to the known form of the 
different ground states, corresponding to the block structure of the A 1 's [TT|[l3] , 
while the second part is the new symmetry-based classification of ground states 
which can be extended to the two-dimensional scenario. 

Proof. First, the G-invariance of the A 1 implies that 

\M{A\U g )) = \M{A\U h U g ul)) = \M(A\U hgh -i)) , 

i.e., all closures from the same conjugacy class are equivalent. Let D 1 , . . . , D J , 
J > I, be all irreducible representations of G. Since the characters Xi(d) = 
tr[D 1 (g)] form an orthonormal set for the space of class functions (i.e. the 
functions which are constant over conjugacy classes), the ground state space 
is equally spanned by 

^triD^g-^MiAlUg)) = \M(A\ ^tr^^" 1 )]^)) ex \M(A\IU)) ; 



Note that IT = for i > I (i.e. if U g does not contain D l ). 

Conversely, linear independence of |.M(A|IT)), i = 1, . . . ,1 can be seen as 
follows: For any class function A g , and with C the L-fold blocking of A, we have 
that 
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Thus, 



= £>|M(^|ni)) = ^M^tr^- 1 )]^!^)) 

i i,g 

i,g i 
= fli V i , 

which proves linear independence. (Note that we had to use that D 1 is contained 
in U g , otherwise IT = 0.) 

The second statement follows from the fact that | Ai(A\U g )) is constant on 
conjugacy classes, and that the number of conjugacy classes equals the number 
of irreducible representations of G. □ 



5 Two dimensions: G— inject ive PEPS 
5.1 Definition and basic properties 

Having understood the one-dimensional case of G-injective MPS, let us now 
turn towards two dimensions. We will introduce some new conventions for the 
diagrams (as in principle, we need a third dimension), which we will explain 
right after the definition of G-injectivity. 

Definition 5.1. Let g i-> U g be a semi-regular representation of a finite group 
G. We call a PEPS tensor A G-injective if 
i) It is invariant under U g on the virtual level, 



(5.1) 



ii) There exists a left inverse to V(A) such that V(A) 1 P(A) = liu, the projec- 
tor on the V g -invariant subspace: 




a? (AT 




(5.2) 



Let us briefly explain the differences in notation: We will try to avoid three- 
dimensional plots as far as possible. PEPS tensors are generally depicted "from 



the top": The four legs in (5.1) are the virtual indices, and a black dot denotes 
the physical index (for the inverse, it is in the lower right instead of the upper 
left corner). As we use the let inverse to make the virtual subspace accessible 
via the physical indices, we will depict only the situation after applying the left 
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inverse whenever possible, as depicted on the very right of (5.2). In order to 



distinguish the "original" virtual level and the one after the application of the 
left-inverse, we shade the latter gray. (This corresponds to the lower and upper 
layer, respectively, in the ID case.) 

Note that there is no need to choose the same representation U g for the 
horizontal and vertical direction. In fact, as mentioned earlier representations 
are assigned to links, and every link can carry its own representation - all that 
matters is that the two tensors acting on a link act with the same representation. 
It is this possibility of changing the representations which enables us to prove 
that G-injectivity is is preserved under blocking. 

Lemma 5.2 (Stability under concatenation). Let A and B be G-injective ten- 
sors. Then, 

•^=^A}^r ( 5 - 3 ) 

(with blocked up, down, and physical indices) is also G-injective with left inverse 



(5.4) 



Proof. First, note that for any two semi- regular representations U g and V g , 
W g = U g (8> V g is again a semi-regular representation. Then, it is clear that C is 
also G-invariant, as the action of the U g on the inner link cancels. 



That (5.4| 



is left-inverse to G follows using tr[U g A] 



Lemma 



4.6 




(5.5) 



□ 



Observation 5.3. Note that G-injectivity is also preserved when contracting 
legs of an already connected block, e.g. the up leg of A with the down leg of B 
in (5.3): The resulting tensor is clearly again G-invariant, and the left-inverse is 



obtained by contracting the corresponding legs of (5.4) with any any operator 



with nonzero trace (e.g., A): The group elements attached to the two legs cancel 
out, as they belong to the same tensor. 



5.2 Parent Hamiltonians 

The idea to construct parent Hamiltonians is essentially the same as in one 
dimension. We define the local Hamiltonian as 1 minus the projector on the 
span of a 2 x 2 block (the smallest block which allows for an overlapping tiling 
of the lattice), and study how the ground state subspace behaves when growing 
the block. For simplicity, we first grow the 2x2 block in one direction until 
we reach a 2 x L lattice, and then in the other direction until we have the full 
L x L lattice with open boundaries. Finally, we study what happens when we 



21 



close the boundaries. While the growing will work essentially exactly as in ID, 
closing the boundaries will give a richer structure. 

Theorem 5.4 (Intersection property) . Let A, B be G-injective. Then, 

1 fJP 



ftEB 



N 



X 



=■■ HM)) 



\P(N)) 



ICW> 



Here, we have chosen to shade the inside of the "boundary condition" tensors 
M, N, and X. 

Proof. The proof is exactly analogous to the one-dimensional case, Theorem |4.8| 
The r.h.s. is contained in the l.h.s., as any can be written as both \a(M)) 

and \(3(N)). Conversely, any element in the intersection can be written as 
\a(M)) — \f3(N)) for some M and N; as in the one-dimensional case, we can 
assume both to be G-invariant. To recover M, we apply the left inverse (5.4 1 
to the left two physical modes of \a(M)) — \/3(N)) and obtain 



j — lJm 





which implies that the state is of the form \((X)). 

Theorem 5.5 (Closure property). For G-injective A, B, C, and D, 



'WL 



HM)) 




m )■ n 



n )■ n 



P )■ n 




"=:|W)> 




□ 



(5.6) 



= ■ 10 



Here, the sum runs over all pairs (g,h) G G x G such that gh = hg. 

Note that if A, B, C, and D arise from blocking the original tensor, this 
corresponds to a closure with two unitaries Uf L and U® L at the horizontal and 
vertical closure, respectively. 
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Proof. The proof again follows closely the proof for one dimension. First, it 
is clear that the r.h.s. in contained in the intersection by choosing M, N, P, 
and Q appropriately. To show that every element in the intersection is of the 
form we consider an element of the intersecion, which can be written as 
\a(M)) = \S(Q)) with boundaries M and Q, respectively. As before, we can 
assume the boundary conditions to be G-invariant. To recover M, we apply the 
left-inverse 



(5.7) 




to \a(M)) = \6(Q)), and obtain 



M 



















a.b.c.d 




= y x, 

£^ h l ,h 2 ,v lr v 2 

h v h 2 ,v v v 2 



v 2 h 1 =h 2 v 1 




(5-8) 

where hi — b a, hi = d 1 c, v± = c a, and vi = d 1 b, which results in the 
constraint t>2hi — h2V\ on the sum. 



Applying the inverse ( |5.7| also to \a(M)) — \f3(N)), we obtain 



i r 

A L 




[note that the left and right block of (5.7) are the inverses for the corresponding 
blocks in |j8(JV))], which proves that hi = hi =: h in (5.8). Finally, the same 



for \a{M)) — \"f(P)) yields v\ — V2 =: g, and the constraint «2^i = h^vi gives 
gh = hg. □ 

Definition 5.6. For a G-injective PEPS tensor A, and g,h € G, define 



\M(A\(g,h))) := 




(5.9) 
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to be the L x L PEPS built from the tensor A, with closures g and h. (The 
representation of g is a property of the bond and omitted for brevity.) 

Theorem 5.7 (Parent Hamiltonian) . Let A be G-injective, 

<?2x2 



a ■ n ■ 



M 



and define 



(i - n 52 



(5.10) 



where the projector Tls 2X 2 ac ^ s on the square {i, i + 1} x {j,j + 1}. Then, 



Hr, 



L 



(5.11) 



has a subspace of frustration free ground states spanned by the PEPS \M.(A\(jg, h))) 
with {g, h)-closed boundaries, where gh = hg. 

Proof. The proof goes exactly as in one dimension: We divide the Hamiltonian 
in four blocks with open boundary conditions (one which does not cover the 
boundary, two which cover the horizontal or vertical boundary, respectively, 
and one which covers both boundaries). For each of the blocks, by virtue of the 



intersection property (Theorem 5.4) the ground state subspace is given by one 
of the sets in (5.6), and thus, Theorem 5.5 yields the desired result. □ 



Definition 5.8 (Pair-conjugacy classes). For (g,h) G G x G, (g',h') e G x G, 
define the equivalence relation 

(g^ti) ~ (g',hf) :<^3xeG: (g, h) = {xg'x' 1 , xtix~ r ) . 

It divides G x G into disjoint equivalence classes 

C[(g,h)] = {(g',ti)\(g,h)~(g',ti)} 

which we call pair-conjugacy classes. To each pair-conjugacy class C, we define 

\M(A\C)):=\M(A\(g,h))) (g,h)eC. 

Note that \M(A\C)) is well-defined, since the c/-invariance of A shows that 
\M{A\{g,h))) - \M(A\(g',h'))) for (g,h) ~ {g',h'). 

Theorem 5.9 (Structure of the ground state subspace). Let C\ U • • • U Ck C 

G x G be the pair-conjucacy classes of G x G for which gh = hg for (g, h) € 
Cfe. Then, the ground state subspace of Theorem \5.T\ is K-fold degenerate and 
spanned by the states \M(A\Ck)) ■ 
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Proof. The proof resembles the one-dimensional case: As observed in Defini- 



tion [OJ closures from the same pair-conjugacy class describe the same state. 
It remains to show that the states corresponding to different pair-conjugacy 
classes are linearly independent. Let 



= J2^\M(A\C k )) =^\ k \M(A\(g k ,h k ))) , 



(5.12) 



with (g k ,hk) G Cfc representatives of C k . Denoting the L x L block of A's by 
G, and applying the left inverse of G, we have that 

k | ^r\H k,x 




with (g k ,h%) € Cfc other representatives of Gfc. This is, (5.12) implies that 
^2^ k Uhi ® Ug* = 0, and since semi-regular representations are linearly inde- 
pendent (Lemma 4.6), this gives X k = Vfc. □ 



It is straightforward to see that the number of different closures, i.e., the 
dimension of the ground state subspace, is equal to the number of particle types 

C[(h,k)], 



A pair-conjugacy class G 



of the quantum double model of G 56 
hk = kh, can be characterized by specifying i) the conjugacy class C[h] of h, 
and ii) the conjugacy class of the normalizer N[h] = {k : hk = kh} of h which 
contains k. Since the number of conjugacy classes of N[h] is equal to the number 
of irreducible representations, this is equal to the number of particle types in 



the G quantum double 57 



Let us now give an intuitive explanation why the closures should be done 
using two closed loops of identical unitaries Uf N and U® N , and why one should 
require gh — hg. The key observation here (which will be essential for the 
introduction of topological excitations later on) is that strings of U g 's can be 
deformed freely due to the G-invariance of A. To this end, consider a string of 
Ug's in the lattice, 





where on the right, we have introduced a new simplified notation which allows 
for more compact diagrams in situations where we care only about the pattern of 
Ug's on the virtual level. Using the G-invariance of the tensors, we can deform 
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the string continuously: e.g., 















i-j 























(5.13) 



and so further. Note that whether g or g 1 has to be used depends on whether 
the (oriented) string of U g 's crosses the (oriented) edges from the right or from 
the left. 

The possibility to arbitrarily deform strings of g's and h's, together with 
the fact that g and h commute, i.e., the that the two strings don't interact, 
implies that we can move the loops anywhere we want on the torus. Thus, the 
closure cannot be detected locally, and the \A4(A\(g,h))) are all ground states 
of the Hamiltonian (5.11). (Note, however, that this does not imply that the 



reduced density operators are the same, but only that they all live within the 
same subspace S2 ■ Indeed, in one dimension the reduced operators can look 
different; e.g., the GHZ state MPS A 1 — \i) (£|, i — 0, 1 gives a Hamiltonian for 
which all a |0, . . . , 0) + /3 |1, . . . , 1) are ground states. 



6 Isometric PEPS 

6.1 Definition an basic properties 

We will now consider a subclass of G-injective PEPS, namely those for which 
V(A) is an isometry (mapping the space of G-invariant tensors unitarily to the 



range of V(A)). Using Observation 2.4 which says that we can w.l.o.g. restrict 



the physical system to the range of V(A), we obtain that a T(A) is isometric if 
V(A)- 1 = V(A*), where the dag ger is with respect to the virtual levels: 

A = J2 \% \ a > & (7. S\ -> ^ := E A UjS I j )p 17. S ) («. P\ ■ 

(Loosely speaking, the reason to use the dagger instead of the conjugation is 
that the edges for V(A) and V(A)~ 1 have opposite orientation.) 

For clarity, we will illustrate proofs for the one-dimensional case where ap- 
propriate. 

Definition 6.1. A G-injective PEPS is called G ' -isometric if U g = L g is the 

left-regular representation, and if V(A)~ 1 = V{A^ ), or equivalently, if V(A) 
restricted to its domain and range is unitary. 

Here, the left-regular representation L g : C' G ' — > C' G ' acts as L g \h) = 
\gh). The following lemma explains why we reqire the restriction to the regular 
representation. 



2G 



Lemma 6.2 (Stability of isometry under concatenation). G-isometry of tensors 
is stable under concatenation. 

Proof. U g = L g A = yjjyl, and thus the left- inverse of the concatenated 
tensor in (4.9) is (B^^A 1 )^ — (A' l B : ')^ - note that the ordering of operators (as 
indicated by the arrows) is reversed for the layer with the inverse P(A)~ 1 . □ 

Note that we do not block indices at this stage, as this would change the 
representation. We will show how this can be done in the section on renormal- 
ization. 



6.2 Equivalence of virtual and physical system 

The importance of G-isometric PEPS lies in the fact that the way the physical 
subspace and the virtual subspace (restricted to the G-invariant subspace!) are 
connected corresponds to a unitary, i.e., a physical operation. This implies that 
any unitary on the virtual level which leaves the G-invariant subspace invariant 
can be implemented by means of a unitary acting on the physical system. 

Lemma 6.3. Let A be a G-isometric PEPS. Then the unitary transformations 
V which can be implemented on the virtual level by unitarily acting on the phys- 
ical system are exactly those which commute with the symmetry, i.e., those for 
which 



(6.1) 





- _ ~<f> 















Proof. Let V be a unitary which satisfies (6.1). Then, V acts unitarily on the 



G-invariant subspace. Thus, the physical operation 





i.e., it acts as V on the virtual indices. Conversely, any unitary operation U on 
the physical level results in a G-invariant unitary V = T'(A)~ 1 UP(A) on the 
virtual level by virtue of 



=:V 
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where we have used (2.4). □ 

A particularly appealing perspective on isometric PEPS is the following. 
Observation 6.4. Every G-isometric PEPS tensor A is isomorphic to 



(6.2) 



by acting unitarily on the physical system]^] Moreover, this property is stable 
under concatenation. E.g., 



AHAh = 



g 



(6.3) 



up to discarding unneeded local degrees of freedom (i.e., restricting the physical 
system to the subspace actually used, cf. Observation 2.4 1. 



Eq. (6.2) follows directly from applying the left inverse (5.2), using that 



it is an isometry. While stability under concatenation can be inferred from 
Lemma 6.2 ("Isometry is stable under concatenation"), it is instructive to see 



how ( 6.3 ) can be obtained from two tensors of the form ( 6.2 ) in a reversible way. 



Concatenating two tensors of the form (6.2), we face a state of the form 







Then, we measure the two sites labeled a and b in the standard basis, and thus 
infer g~ 1 h = ab^ 1 . (We will use a b, etc. for both the site and its state.) Now, 
apply U g-^h to sites c, d, and e. This maps Uh i-> UhW_ lh = U g (on site c) or 
f//i-i <— > Ug-ihUh-i —Ug-i (on sites d and e, as they are bras), and thus yields 



a tensor of the form (6.3) 



Do we loose reversibility in this process, given that we measure? No: One 
way to see this is to replace the measurement followed by Ug-i^ by a unitary 
controlled by the state of a and b, which keeps the process reversible. In par- 
ticular, the state of a and b is separable from the rest after this operation (it is 
the only state which still depends on h). This, however, implies that we could 
equally well discard (i.e., measure) it, and rebuild it if we need to reverse the 
mapping. 



4 Note that this decomposes the physical level into four systems, represented by kets or 
bras depending on the directions of the arrows. 



28 



6.3 Blocking and renormalization 

An aspect which has to be taken care of separately for G-isometric PEPS is 
blocking. Recall that the proofs for the 2D case involved taking a block of G— 
injective tensors A with representation U g and regarding it as a new G-injective 
4-index tensor B with representation V g — U g <S> • • • <£> U g . This was possible 
since a tensor product of semi-regular representations is again semi-regular. In 
contrast, the regular case needs some special attention. 

The main observation to deal with blocking in the regular case is to observe 
that for the left-regular representation L g , 

(L g )® N ® (L\) 9M = L g ® . 

(This follows since both representations have the same character, which means 
they are isomorphic.) This implies that any number of bonds with a regular 
representation can be mapped - on the virtual level - to only one bond carrying 
the L ff -symmetry, whereas all other bonds have no symmetry. 

Let us now first show that this isomorphism can be implemented by a physi- 
cal operation, before discussing what this implies for the structure of the PEPS. 
For simplicity, we restrict to the case of two bonds, L g ® L g = L g ® 1. The 
isomorphism is implemented by the map 

T = ^2\a,ab){a,b\ : T ] (L g <g> L g )T = L g ® 1 . 

a . 6 

Let us now consider two tensors A and B adjacent to the L g ® L g bond. By 
virtue of Observation |6.4[ and grouping all other indices, they are isomorphic 
to 





_ a_ 


*® 






h,k 




*® 







Now, we apply the unitary T (T^) to the pair of indices labelled a, b (c, d): This 
physical transformation changes the network to tensors A', B' with symmetry 
L g (g) 1: 

Kk J ■ J /a Arf _ Bd " 

As illustrated on the r.h.s., this means that A' and B' are actually a tensor 
product of two independent tensors: A = A u <E) A& and B = B u ® B^. Whereas 
the upper tensors have an L g -symmetry, the lower tensors have no symmetry, 
and form a maximally entangled pair between adjacent sites. Thus, these de- 
grees of freedom do not contribute to the non-local properties of the PEPS, and 
can thus be discarded. 

Observation 6.5 (Blocking indices preserves G-isometry). Blocking indices of 
G-isometric tensors gives a new G-isometric tensor, up to local isomorphisms 
and discarding maximally entangled pairs between adjacent sites. 
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Let us now show that this procedure establishes a renormalization scheme for 
PEPS. The idea of renormalization is to study non-local properties of quantum 
many-body states by grouping several sites into one, identifying and discarding 
"irrelevant" (i.e., local) degrees of freedom, and iterating until convergence is 
reached. 

Consider now a G-isometric PEPS with tensor A (and symmetry L g ), and 
consider a block of 2 x 2 tensors. The blocked tensor has then symmetry L g (g)L g , 
which by virtue of the above procedure can be mapped to a L g ® 1 symmetry 
by local operations: 




Here, the / tensors are isometric tensors with no symmetry, i.e., they map 
the virtual system one-to-one to the physical system. Thus, by local unitaries 
on 2 x 2 blocks, the G-invariant PEPS with tensor A has been transformed 
into a tensor product of the following states: First, the same PEPS described 
by the tensor A, but on a coarse-grained lattice, and second, the state given 
by the identity tensors I - a tensor product of maximally entangled states 
between neighboring blocks, which can be discarded if only looking for non- 
local properties. This results in the following observation. 

Observation 6.6. G-isometric PEPS are renormalization group fixed points^] 

A similar procedure should also apply to the case of isometric tensors with a 
non-regular representation U g of the symmetry, or even to non-isometric PEPS: 
(Approximate) decompositions of the tensorized symmetry Uf k will yield (ap- 
proximate) RG schemes for the PEPS. Another application of this RG scheme 
will be illustrated in the examples section: It can be used to get rid of extra 
symmetries in the original tensors which are of local nature and vanish after a 
step of blocking. 



6.4 Structure of the ground state subspace 



We will now study the properties of PEPS with a (g, /i)-closure, Eq. (5.9), which 
appear as ground states of the 2D parent Hamiltonian (5.11), for the case of 



G-isometric PEPS. We want to understand three things: First, is it possible to 
map between arbitrary states in the ground state subspace by acting only on 
a restricted region? Second, are different ground states locally distinguishable? 
Third, what is the entropy of a continguous block of spins? 

All these questions can be addressed by considering the following scenario 



5 Since the regular representation is the only faithful one verifying U g ®U g = U g <&i, G— 
isometric PEPS seem to be the unique fixed points of this type of renormalization procedure. 
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with G-isometric B and D: 




(6.4) 



which by observation 6.4 is locally (where local refers to the regions B and D) 
isomorphic to 

(h) 



(6.5) 




When dividing a large lattice in two blocks as above, we need to include an RG 
step in which we separate out 3D — 4 maximally entangled states between B 
and D; here, the boundary 3D is the number of bonds crossing the boundary 
of D before the RG step. 

We will now devise a unitary transformation acting only on the B region, 
which will decouple h and g from D. To this end, consider the sites marked 
a, 6, and c above. Since all group transformations are left regular, we have 
c = x~ 1 gxb ) and thus cb^ 1 = x~ 1 gx: Having access to b and c actually gives us 
access to x~ 1 gx. Thus, a unitary transform 

|c) (6| (a\ ^ |c) (b\ (cb- 



g xa 



)\b) 



(6.6) 



will remove the left g in (6.5) (the one marked *) from the boundary. A corre- 
sponding transformation allows us to remove the lower h, leaving us with 



x,y 





® , (6-7) 



where in the second step we have substituted z := x^ 1 y. Now, however, g, h, 
and x only act on degrees of freedom which are completely inside the region 
B, while the degrees of freedom of D are equal to the corresponding ones of B, 
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up to a constant shift z. In formulas, the state above can be written (up to 
isomorphism) as 

®ICG7,ft)> B , 

z 

with 

H Z ))BD = J2\ zb )B \ b )D > 
b 

\C(g,h)) B = \x~ l gxp) B \p) B \x~ 1 hxq) B \q) B . 

x,p,q 

Note that additionally, there are dD — 4 maximally entangled states between B 
and D which have been separated in the RG step. 



Equation (6.7) has several immediate consequences. 



Theorem 6.7 (Local equivalence of ground states). For G -isometric PEPS, 
it is possible to unitarily transform between any two states in the ground state 



subspace ^2 C Ac\-M.(M\C)} of Theorem 5.9 by acting only on two stripes (each 
of width one) wrapping around the torus. 



Proof. By fixing two such stripes, the PEPS is partitioned as in (6.4), with the 
stripes labelled B and the (topologically trivial) rest D. Then, decomposition 
(6.7) shows that the states \A4(M\(g,h))) spanning the ground state manifold 
differ only within B, allowing to unitarily transform between any two states by 
acting only on B. □ 

Corollary 6.8 (Local indistinguishability of ground states). For G -isometric 
PEPS, the states in the ground state subspace of Theorem \5.9\ cannot be distin- 
guished by local operations, i.e., those which act only on a topologically trivial 
region. 

Proof. This follows immediately by using that for any topologically trivial re- 
gion D, it is possible to transform between any two states in the ground state 
subspace by acting with a unitary on its complement, B. □ 

Theorem 6.9 (Topological entropy of G-isometric states). For a G-isometric 
PEPS, the reduced density operator pjj of any topologically trivial block D has 
rank |G| 9D_1 and a flat spectrum. Here, the length dD of the boundary of D 
is defined as the number of bonds crossing the boundary. Thus, both the von 
Neumann entropy S and all Renyi entropies S a of po are 

S( PD ) = S a ( PD )=log\G\(dD)-\og\G\ , 

where log \G\ is the topological correction to the area law l5S\\6(Jj. 



Proof. We can again use the partitioning (6.4): We can then infer from (6.7) 
that (after unitary transformations on B and D, and discarding the unentanglcd 
part of B), the state is of the form |/z(z))® 4 , with 

\^(z))^J2\ zb )B\b) D ■ 
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Additionally, we also have (dD — 4) maximally entangled states from the 

initial RG step. It is now possible to use one copy of \p>(z)) as a reference frame 
which allows to remove the reference to z from the other copies: Take two copies 
of and apply the operation \p, q) i— > \p,p~ 1 q) both on B and G. Then, 

\fx(z)) \fx(z)) = ]T \ zb > zc )b % c) D b ~ 1c )b l 6 < b ^ c ) D = IM*)) • 

be be 

Thus, all but one \ti(z)) can be mapped locally to a maximally entangled state 
of Schmidt rank |G|, while the remaining "reference frame" copy - the only 
z-dependent contribution to the state - turns into ^ z which is separable. 

Thus, B and D share a total of {dD — 1) maximally entangled states, which 
yields a flat Schmidt spectrum of rank |G| 9D ~ 1 □ 

Corollary 6.10 (Topological Renyi entropy of G-injective states). Let A be G- 
injective, with U g = L g the regular representation. Then, any reduced density 
operator pu of any topologically trivial block D has rank \G\ dD ~ 1 , and thus its 
zero Renyi entropy is 

S { PD ) = \og\G\{dD)-\og\G\ . 

Proof. This follows immediately from the preceding theorem, since any G-injec- 
tive A with U g = L g can be transformed to a G-isometric one by a local (non- 
unitary) transformation on the physical system. (This can be seen by doing 
a singular value decomposition of V(A) — UDV: UD is the desired linear 



operation, which by Observation 2.4 has full rank.) □ 



6.5 Commuting parent Hamiltonians 

Let us now show that the parent Hamiltonians of G-isometric PEPS are special: 
They are sums of commuting local terms. We will demonstrate the proof for 
MPS, but the generalization to PEPS is straightforward. 




Lemma 6.11. For G isometric MPS, the local terms hi [Eq. (5.10)/ in the par- 
ent Hamiltonian - which project onto the complement 0/S2 — tr[A*A 3 'X] \ij) x\ 
- are of the form 

(6.8) 

6 Note that we could have used the same argument without the initial RG step: In that 
case, we would have arrived at dD copies of \fi(z)) which could have been converted to (dD — 1) 
maximally entangled states the same way. 



33 



Proof. First, hi is a projector, since 




it J 



Second, its range is clearly contained in 1S2, and third, any state of the form 
J2ij ti^A'X] is preserved by he 





A 




A 




r 


[ 


A" 




A f 


: 














r 


A 




A 








X 








X 



Finally, hi is self-adjoint, which proves the lemma. 



□ 



Note that the corresponding operator (6.8) defined for non-isometric PEPS 
(with A^ 1 instead of A^) fails the last condition: while it is a projector onto 1S2, 
it is not self-adjoint, i.e., not a Hamiltonian. 

Theorem 6.12 (Commuting Parent Hamiltonians) . For G-isometric PEPS, 



the terms hi of the parent Hamiltonian of Theorem \5. 7| cf. (6.8), commute. 
Proof. We start from the identity 




(Note that this the projector onto the joint kernel of hi and /12, which is equal 
hih2 for commuting hi, /12O The l.h.s. can be transformed to 



g 



A Kg)] A 



g 





i.e., hih2, while the corresponding transformation on the r.h.s. yields /12/11, 
proving that /11/12 = ^2^1- The same proof applies to two dimensions. □ 
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6.6 Anyons 



In the following, we show how to understand anyonic excitations of G-isometric 
PEPS using the symmetry of the tensors. Intuitively, excitations should be 
formed by (open) strings of U g 's, similar to the (closed) strings distinguishing 
the different ground states: Since such strings can be continuously deformed, 
they cannot be detectable anywhere locally, except for the endpoints. As we 
will show, these strings give rise to one type of particles: fluxes, but there will 
also be a second type of complementary particles: charges. 

Note that in the following, we consider the bulk of a PEPS formed by G— 
isometric tensors A, and assume the boundaries to be far away: The results 
generally hold independent of the form of the boundaries. 



I. Magnetic fluxes 

Definition 6.13 (Fluxons). Consider a PEPS formed of G-isometric tensors 
A, let g G G, and let U g = L g be the left-regular representation. A pair of 
magnetic fluxes or fluxons in the state (g, g~ l ) is defined by placing a string of 
U g and U g -i on the virtual level as follows: 



(6.9) 



The fluxons are attached to the two marked plaquettes forming the endpoints 
of the string; they are characterized by the fact that there is an odd number of 
g and g^ 1 on the adjacent edges. Whether g or g^ 1 is used depends on the 
orientation of the string relative to the bonds; we will generally talk of a "string 
of g 's" meaning both g and g^ 1 . Also, we assign the state g to the starting point 
and g~ x to the endpoint. The particle type of the two fluxons is given by the 
conjugacy classes C[g] and G[g _1 ], while the exact states g and g~ l determine 
the internal state of the fluxons. We will show in the following why these choices 
make sense. 

Lemma 6.14 (Deformations of the string). The endpoint plaquettes of a string 
are fixed. Except for that, the string connecting a pair of fluxons can be deformed 
at will without acting on the system. This implies that it cannot be observed 
anywhere except for the endpoints. 



Proof. The string can be defomed using the G-invariance of the tensors, cf. ( 5.13 1 
as long as the endpoints are not involved. On the other hand, the G-invariance 
cannot change the parity of g's around a plaquette, which implies that the end- 
points cannot be moved. (However, the string can be deformed to reach the 
endpoint from any side.) Since according to the following theorem, endpoints 
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can be measured, there also cannot be any other way to move them away without 
acting on the physical system. □ 

Theorem 6.15 (Detection of fluxons). The particle type C[g] of a fluxon can 
be detected by a measurement on the plaquette supporting it. 

Proof. As we can arbitrarily deform the string of g's ending at a given plaquette, 
it suffices to consider the following situation: 



(6.10) 



We will use the four sites of the plaquette to build a "tweezer" which allows 
us to "grab" g, up to a twirl. To this end, block the four sites to a U-shaped 
tensor. From Observation |6.4[ we know that there is a reversible operation 
which transforms the blocked tensor to 





(6.11) 



By measuring ab" 1 — xgx~ l , we can now infer the particle type C[g] of the 
fluxon. Since the value of x cannot be determined, this is all the information 
we can gain about a single fluxon. 



Looking at (6.11 1 and the way our measurement of C[g] works, one can see 
that we actually do not use any of the outgoing indices, i.e., those which connect 
to the virtual level and thus to neighboring tensors. Therefore, we can simplify 
notation by omitting these indices together with the sum over x and analyze 
the action of the detection tweezer on 



where x is unknown. We will make use of this simplified notation from now on 
in the analysis of subsequent constructions for detecting, moving, and creating 
anyons. 

Note that the reason we built a tweezer and not just acted on the site left 



and right of g in (6.10) is that we needed to synchronize the twirl; otherwise, 
we could have only determined xgy -1 , which is completely random. Let us also 
note that the tweezer can open in any direction, as the string of <?'s can be 
attached to any egde of the plaquette. □ 
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Theorem 6.16 (Moving of fluxons). Fluxons can be moved by local unitaries, 
without knowledge of their state. 

Proof. We give the construction for the scenario where we want to move a string 
to the right, i.e., achieve the following transformation 




by means of local unitaries. To this end, we will again build a tweezer (opening 
towards the left) which will allow us to access the two vertical edges on the 
left. Again, using Observation |6.4[ and neglecting unneeded bonds, the above 
transformation is equivalent to the transformation 





for arbitrary x. This can be accomplished using the same transformation already 
used in (6.6 1 for this purpose: '" 



c i — ^ 



|foc- 1 a> \b) 



(Note that each of the 
kets and bras describes the state of a separate quantum system: \c) (a\ (b\ thus 
describes a pure state of the three systems a, b, and c. Whether we use kets or 
bras is determined by the arrows associated to the links.) □ 

Theorem 6.17 (Creation of fluxons). For any g G G, the fluxon pair 



^2\{zgz 1 ,zg x z 1 )) 



(6.12) 



can be created deterministically by local operations^ Note there is one such state 
per conjugacy class C[g\. 

Proof. We want to implement the transformation 




7 The reason that we can only create the equal weight superposition is that all other 
superpositions carry non-zero charge, which is a conserved quantity, cf. [571. 
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To this end, we build a tweezer on the lower plaquette. Neglecting unneeded 
degrees of freedom, the transformation we want to achieve is equivalent to 




(6.13) 



for some unknown x. This can be accomplished by replacing the state \r) (r\ b 
with ^ s z \zgz~ 1 (s\ b , which can be done unitarily if wanted. □ 

II. Braiding of fiuxons 

In the following, we investigate what happens to fiuxons when they are braided 
around each other. The general scenario is as follows: Consider two pairs of 
fiuxons, and move one fluxon of one pair around one of the other one: 




In principle, it does not make sense to ask how the anyon B\ is affected when 
it crosses the line between A\ and A 2 , since this line can be deformed at will. 
However, it does make sense to think like that if the pair B\ and B 2 is brought 
together after braiding to measure the joint flux, since then Bi eventually has 
to cross this line. We will discuss how to measure the joint flux later on. 

Instead of investigating what happens when we move the fluxon B\ across 
the string connecting A\ and A 2 , we rather keep B\ fixed and move the string 
across it - this is easier to analyze, since it takes place purely on the virtual 
level: 




(6.14) 



Let us now consider this scenario - moving the A\-A 2 string accross the anyon 
B\ - in the PEPS representation on the virtual level. To this end, let A\ and 
A 2 be in the state (ff,g -1 ), and B\ and B 2 in state (h, h^ 1 ). Then, using the 
G-invariance of the tensors, 
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i.e., the Bi end of the string is conjugated by g, and is now in the state ghg~ l . 
Note that this does not change the conjugacy class, i.e., the particle type. This 
is crucial since the crossing of B\ and the A\-A 2 line cannot be assigned a well- 
defined position, and thus the change of the state of B 2 cannot be detected by 
measuring B 2 alone. 

However, when we bring B\ and B 2 together, it is possible to measure the 
particle type Clghg^ 1 ^ 1 ] of their joint flux, and thus infer information about 
the state g of A\. For two fluxes in states a and b, their joint state is given by 
ab; this is consistent with how the two fluxes together affect another flux when 
braiding. Note that actually merging the two fluxes will require a measurement, 
since this is an irreversible process. 

Let us now describe how to measure the joint flux. To this end, consider 
two fluxes a and b sitting on adjacent plaquettes - this is as close as we can 
bring them using the "fluxon moving tweezer" of Theorem |6.16| - and route the 
virtual strings as follows: 




Now, we can use basically the same tweezer construction used for measuring 



single fluxons in Theorem 6.15 to measure the C[a6]: To this end, we build a 
tweezer covering both plaquettes, and connected at the left side, and measure 
C[ab] across the right vertex. In fact, this construction can be extended to 
measure the joint flux of any number of anyons enclosed in the tweezer. 



Note that the braiding (6.14) does not only change the state of Bi, but also 
that of Ai if it brought back to interfere with A 2 . The reason is that even after 
repairing the B particles, the loop of /i's remains there and affects any anyon 
crossing it. Note that while the way A\ is affected when moving it out of the 
loop depends on whether it is moved above or below the path it came from, this 
difference cannot be observed when repairing it with A 2 , and vanishes when A\ 
and A 2 are aligned in a fixed way. 



III. Electric charges 

Let us now define electric charges. Different from fluxes which always come in 
pairs, it is possible to define (though not create) electric monopoles. 

Definition 6.18 (Electric charges). An electric charge ("chargeon") with charge 
c, where c labels an irreducible representation D c of G, is given by a deject 




(6.15) 
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which is attached to an edge of the lattice, i.e., instead of the identity, the 
edge acts as J2 g XgiPd) IsXffl ■ Here, Xc(') = tr[D c (-)] is the character of the 
irreducible representation c, and p 6 G characterizes the internal state of the 
chargeon. 

Theorem 6.19 (Detection of chargeons). The charge c of a chargeon can be 
detected by measuring across the two vertices adjacent to the edge supporting 
the charge. 

Note that we will later show that the internal degree of freedom p is changed 
by strings connecting fluxons, which rules out that we can measure p itself. 

Proof. We will perform a joint measurement on the vertices left and right of 



the chargeon in (6.15). From Observation 6.4 the state of these two sites is 



isomorphic (up to unneeded links) to 



(6.16) 



i.e., we have access to the state 



^2 *c{pg) \xg) (yg\ = ^ Xc{px 1 k) \k) (zk\ 



for some x and y, with k := xg and z 



yx 



As a consequence of the group 



orthogonality relations, = J2k Xc(hk) \k) are orthogonal vectors for different 
irreducible representations c, (v^v^,) = f(h, h')6 CtC i. Therefore, the charge c (as 
well as z) can be determined by measuring the two sites in (6.16). Note that on 
the other hand, no information about p can be learned, since x is unkown. □ 



Theorem 6.20 (Moving of chargeons). 
ations. 



Chargeons can be moved by local oper- 



Proof. Since chargeons are localized perturbations on the virtual level (unlike 
fluxons which arise from strings), moving then is as simple as swapping: In 
order to implement the transformation 




(6.17) 



we will need to block the left two and right two spins independently. The 
transformation is then locally equivalent (up to unneeded degrees of freedom) 
to 

I® ©I |©>W@| ' 
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and thus, by swapping a with b and c with d, the chargeon is moved as in (6.17) 



Note that a similar construction can also be used to move a fluxon around the 
corner. □ 

Theorem 6.21 (Creation of chargeons). For any irreducible representation c 
and any p € G, the chargeon pair described by 



n c , P = Xciph^g) Is) (g\ u <8 \h) (h\ d 

g,h 



(6.18) 



can be created by local operations. 

Note that this gives a particle with charge c in position u and a particle of 
charge c in position d, where c labels the complex conjugate of representation 
c. 



Proof. First, observe that the pair (6.18) is invariant under conjugation by x € 
G: 



(U x ® £7 x )n c>p ([/t ® Ul) = Xciph x g) \xg) (xg\ ® |a:ft) (xh\ 

= Xc{ph- l xx- l g) \g) (g\ ® = n c , p 



(6.19) 



To create the state (6.18), we need to implement the transformation 




To this end, we use a U-shaped tweezer encompassing all six sites, which leaves 
us with the task of implementing the transformation 






(6.20) 



for some unkown x. Due to the [/x-invariance of II C! p [Eq. (6.19)], neither 
side depends on x, and the transformation can be implemented (unitarily) by 
mapping the initial state of the four sites a,b,c,d to H c . p . □ 
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IV. Braiding of charges with fluxes 

Let us finally see what happens when we braid chargeons with other particles. 
Since chargeons are localized defects in the lattice, they do not affect other 
particles which are braided around them. On the other hand, they are affected 
if braided around a fluxon. Consider again the braiding (6.14), now with A\-A 2 
a pair of fluxons in state (k, fc -1 ), and B\-Bi a pair of chargeons with charges 
c and c, e.g. as in (6.18). Again, we will choose to move the string connecting 
the fluxons A\ and A2 over the chargeon B\, rather than moving B\ across the 
string. Analyzing the situation for a single projection \g) (g\, we find that 





This is, \g) (g\ ^ U k \g) (g\ u£ = \kg) (kg\, and thus 

Xc(pg) \g) (s| >-> X! Xc(pk~ l g) \g) 



but 



braiding a chargeon around a fluxon of flux k leaves its type c invariant, 
multiplies the parameter p characterizing its internal state by 

To interferometrically detect this change, one can e.g. start from a pair of 
charges in state H cp (6.18), braid them around the flux k, and bring them back 
together. Then one can reason as follows. We have seen that, by acting with a 
unitary W on six sites, the original state (before the braiding) decomposes as 

^ c >p)abcd 



. \ip x ) where \G\ |U 



c ,PUbcd 



is exactly the right hand side of (6.20) 



which is indeed independent of x. The group orthogonality relations assure that 
|n c .p) abc(J is normalized, which in turn gives Y2 X y{ipx\tpy) — !• Now, after braid- 



ing and applying W we are left with J2 a 



IT 



c,p,k 



1 bed 



\if) x ), where \G\ 



n 



c,p,k 



is the right hand side of (6.20) but with the coefficients given by Xc{ph 1 k 
Then, we can implement in abed (after having applied W) the measurement 



abed 



given by |n c , p ) (II, 



and its orthogonal complement, which outputs that 



the state is unchanged with probability 

Xc{k) 



where \c\ denotes the dimension of the representation. To see that, it is enough 
to use the group orthogonality relations once more. 



7 Examples 

In the following, we present examples of G-injective PEPS. We start by dis- 
cussing Kitaev's code state 56 , including an explanation of how to construct 
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Figure 1: The lattice for Kitaev's code state and the double models. 



the original PEPS and to obtain the G-isometric PEPS by an RG step. We 

, the 
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then consider the generalization of the model to arbitrary finite groups 
so-called double models, and show that they correspond to the set of G-isometric 
PEPS. Finally, we show that by using semi-regular representations, it is possible 
to construct PEPS which are locally equivalent to the double models, yet can 
be realized with a lower bond dimension. 



7.1 Kitaev's code state 

Kitaev's code state can be defined on a square lattice with qubits {|0) , |1)} 
attached to its vertices, and two types of plaquettes, A and B, see Fig. [I] (The 
arrows will be used in the non-abelian case). The Hamiltonian consists of local 
terms associated to plaquettes, each acting on the four surrounding qubits: For 
each A plaquette pa , define 

h p / := |(1 - Z® 4 ) (7.1) 

acting on the four sites adjacent to the plaquette pa', this is, h p A A is the projector 
whose ground state subspace is formed by the configurations \g, h, k, 1} with even 
parity around the plaquette. For each B plaquette pb, define 



= \g + s,h + s,k + s,l + s) (g, h, k, l\ , 



(7.2) 



where addition is modulo 2; this is, h p B B is a projector whose ground state 
subspace is invariant under flipping all spins adjacent to ps- Now let 

h a ■= h A A > and h b ■= Yl h B B ■ 

PA Pb 



Then, Kitaev's code state Hamiltonian 56 is given by H = Ha + Hb- 



Let us now explain how to construct a PEPS representation of a ground 
state of H 44 . To this end, observe that all terms h? A A , hF^ in the Hamiltonian 



commute and are projections: Thus, the ground state subspace is given by the 
product 

nia-^IIK 1 -^)' 
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and a product of local operators can always be written as a PEPS. 

A particularly appealing way to construct the PEPS representation is the 
following. We start from the state |0, ...,0), which has even parity around 
each plaquette and is thus a ground state of Ha- By sequentially applying all 
projections |(1 — h?g) to the initial |0, . . . , 0) state, we thus end up with a 
ground state of HPj 

The projector |(1 — h P g) can be expressed as a tensor network of the form 



Pout 




which takes a layer of physical inputs p m at the bottom, and maps them to a 
new layer of physical outputs p out at the top. The four PEPS tensors above are 
all equal, 

a . k 



where v denotes the virtual indices. The natural way to think of this construc- 
tion is to assign the value s to the plaquette rather than to the bonds: When 
contracting the virtual indices around the loop, there only remains a single sum 
over s, which implements the projection |(1 — h P g) by simulateously flipping 
all four physical spins conditionally on s. 

The B plaquette projections can be divided into two disjoint layers, the 



PEPS representation of each of which is given by (7.3 1. When combining the 



layers, at each site two P tensors (7.3) are stacked atop of each other, with the 



virtual indices pointing in opposite directions. In order to obtain the PEPS 
description of the ground state itself, we have to apply these two layers of P to 
the initial state |0). This yields the tensor (cf. [44]) 



r.s 



T = > '\r + s) p \r,s) v (r,s\ v 

[|0,0)„(0,0|„ + |1,1)„(1,1|„ 



1 



|0,1>„<0,1|„ + 11,0)^1,01, 



(7.4) 



Here, the r's and s's on the virtual level are placed such that they are associated 
with one B plaquette each, 

aIb 

' 

B JA 

and correspondingly rotated for the other type of vertices. 

It is again convenient to think of the virtual indices as assigned to B plaque- 
ttes, rather than to bonds. Then, the value of each qubit in the lattice is given 



8 Note we have to make sure that |0, . . . ,0} has non-zero overlap with the ground state 
subspace of Hb ■ This is true as it has non-zero overlap with (|0) + 11)) (§> • • • ® (|0) + |1}) which 
is in the kernel of Hg. 
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by the difference r — s = r + s between the virtual degrees of freedom r and s 
assigned to the two adjacent B plaquettes. Alternatively, one can think of the 
state as arising from |0, . . . ,0) - an eigenstate of Ha - which is symmetrized 
by coherently flipping the spins adjacent to any B plaquette, controlled by its 
virtual degree of freedom. 



From (7.4 1, it is straightforward to see that T is characterized by the follow- 



ing symmetries: 




Thus, T has the desired Z2-invariance with respect to X® 4 , but it is lacking 
Z2-injectivity since it has additional symmetries with respect to Z® 2 . In order 
to get rid of these symmetries, we consider a 2 x 2 block around an A plaquette, 
which has symmetries 




The X® s symmetry is the topologically interesting one, whereas the Z® 2 sym- 
metries, after blocking, only act between adjacent tensors: They should there- 
fore be local in nature. Moreover, since X ® X and Z ® Z commute, there has 
to be an isometry which separates their action: This is accomplished by the 
CNOT gate \a) \b) i-> \a) \a+ b) which we have already used for the RG flow in 
Sec. |6.3[ and which maps 



X ® X i-> X ® 1 
Z®Z^ 1® z . 

By applying these CNOTs, the 2 x 2-block of T's is thus transformed tc^] 

ft 





where the tensors K and L have exactly the symmetries 




and 



9 Note that since V(T) is isometric on its support, the RG step can be implemented by a 
physical unitary and a relabelling of bonds. 
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respectively. Thus, the tensor K is the desired Z 2 -isometric tensor describing 
the topological code state, while the L are tensors with only one virtual qubit 
index with Z2 symmetry: They therefore give rise to product states between 
adjacent sites. Note that this also explains why the original PEPS representation 
of the code state needed double as many bonds as actually required to obtain 
the topological entropy: The original tensor T has extra local symmetries. 

Let us now try to see how the blocked and renormalized tensors K actually 
look like. As compared to the original lattice, the tensors sit on top of (i.e., 
contain) an A plaquette, and each bond corresponds to a B plaquette. The 
other half of the A plaquettes forms the plaquettes between the renormalized 
tensors: 




The summation index s in the PEPS construction (7.3 1 and (7.4), which we 



argued was associated to B plaquettes, is now associated to edges, i.e., it became 
an actual bond state. Going through the RG scheme (or observing that the 
action of the tensor P, Eq. (7.3), is to add s to all qubits surrounding the 



associated B plaquette), we find that 



K = ^2 \p + q, q + r, r + s, s + p) p <g> \p, q) v (r, s\ v 
\1 



pqrs 



E p 

pqrs 



(p+q) 



(7.5) 



One possible perspective on the PEPS representation is thus to assign values 
or 1 ("colors") to each B plaquette, and let the physical states be the difference 
between the colors of the adjacent plaquettes; the state is then given by the 
superposition over all color patterns. Note that V{K) does not have full rank, 
as the physical subspace in (7.5) is restricted to states with even parity. 

What is the parent Hamiltonian for the blocked tensor? As it turns out, it 
can be split into three components, each of which can be associated to a term in 
the original Hamiltonian: i) a local term which enforces the local constraint; ii) 
a four-site term enforcing the plaquette constraint for the A plaquette covering 
four tensors; and iii) a two-site term making sure that the state is an equal 
weight superposition of all edge colorings: it enforces that the state is invariant 
under adding 1 to all four qubits surrounding the edge. 
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As we have seen earlier, there is also a different way to express the PEPS 
tensor K, using Observation 6.4 




Written in this form, the four physical qubits are +1 eigenstates of X® 4 , rather 



than Z® 4 as in (7.5), and the different terms in the Hamiltonian ensure i) that 
this local constraint is obeyed, ii) that the product of twirls around a plaquette 
vanishes, and iii) that all possible twirls on an edge appear with equal weight. 



7.2 The double models 

Kitaev's code state Hamiltonian can be naturally generalized to any finite group 
G. To this end, identify the states at each site with group elements \g), and 
generalize the Hamiltonian terms as follows: 

h A :=1 - E \9,h,k,l){g,h,k,l\ , (7.6) 

ghkl— 1 

where the sum is over all group elements whose product ghkl equals the identity 
1, and the product is taken in the order indicated by the arrows in Fig. [I] from 
left to right; and 

KP B : = 1 - }Z \ U >9> U ' h ' Usk > U ^ & h ' fc ' 'I ' ( 7J ) 

s 

where U s = L s (U„ = R s ) if in Fig. [l] the arrow on the spin points away from 
(towards) pb', we keep this convention throughout. Here, L s : \s) n> \sg) and 
R s : | s) i — >- |<7S -1 ) are the left and right multiplication, respectively. The total 
Hamiltonian is again the sum of these terms for all plaquettes. 

The PEPS construction can be carried out the same way as for the Z2 case, 
by starting from |1, . . . , 1) and applying projectors 

s,k 

applying two of these projectors to |1), we arrive at the PEPS tensor 

T = Y / U r U s \l) p \r,s) v (r,s\ v . (7.9) 

r,s 

Here, the action of U r and U s is determined by the site it acts on, and by the 
plaquette associated with r and s, respectively (which implies the order does 
not matter); U r U s |1) is thus either |rs _1 ) or |sr _1 ). 
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The tensor T has again additional symmetries, which can be removed by 



blocking around an A plaquctt 



this yields the G-isometric tensorfor the quantum double, 



and renormalizing via |o) \b) i— > \a)\b 1 a); 



A = \ pq 1 ' qr 1 ' rs 1 ' sp X ) P ® \p>i)v ( r ' s L 



z ^ 

pqrs 



(7.10) 



Note that this PEPS representation can again be interpreted by assigning colors, 
i.e., group elements, to the bonds, and letting the state of each qu-|G|-it be the 
difference between the adjacent bonds. 

The Hamiltonian acts analogously to the Z2 case: i) a local term ensures 



that the four sites in (7.101 multiply to the identity; ii) a plaquette term ensures 
the same for the corresponding sites of four tensors, and iii) a two-body term 
enforces equal weight for all bond configurations by applying U s coherently to 
the four adjacent sites. 

Alternatively, we can again use the representation 




= z<zy-®~ , (7.11) 

S 



where the role of the Hamiltonian terms in this representation are to ensure 
that i) all £/ s -invariant local configurations have equal weight, ii) the product 
of the g's around a plaquette is the identity, and iii) for any bond, all twirls 
appear with equal weight. 



7.3 Reducing the bond dimension: semi-regular represen- 
tations 

In the following, we would like to go beyond regular representations, and show 
how a PEPS which is G-injective with respect to a semi-regular representation 
can implement topological states which are equivalent to the quantum doubles, 
but with a lower bond dimension. 

To this end, consider the tensor K in the representation ( |7.11| ). Here, g is 
the left-regular representation of a group, and can thus - by a proper relabelling 

10 The PEPS representation depends on the rotation direction prescribed by the arrows 
around the plaquette. Eq. ||7.10| is for clockwise rotation. 
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of the virtual indices - be chosen to be of the form 



Here, D l , i = 1, . . . , I, are the irreducible representations of g, with dimensions 
di. Let us now take the smallest semi-regular representation 



V a =($D\g) 



of G, which contains every irreducible representation with multiplicity one, and 
define 

i=i 



[Note that 6 4 = A, with A from (4.6 1.] With this, we define a tensor 



Sh=L 




(7.12) 



We now claim that the state given by the tensor H is the same as the double 
model for the group G, i.e., the state described by (7.11) with the regular 



representation, up to unitaries between neighboring tensors which act on disjoint 
subsystems, i.e., local operations. This demonstrates that the double model has 
unneeded local degrees of freedom, and the bond dimension needed to represent 
the model is J^i <k rather than J^, d\. Note that this is different from the bond 
dimension gained by renormalization in that no blocking of tensors is required. 

In order to prove local equivalence of the two states, we will explicitly con- 
struct an isometric transformation which maps between the two states. (It has 
to be isometric as the dimension of the underlying space changes.) The trans- 
formation will act on the two endpoints of a bond simulateously, where it needs 
to implement a mapping between 



and 



this is, between 
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and 



X! ^2 D ki,k(9 h ^Hikim) {i;h;mi\ 



i—l ki,L—l rrii — 1 



This, however, can be accomplished by mapping between the basis elements of 
the two spaces as 



\i;ki) (i;h \ ^ 



-= ^ (i;k;mi\ ; 



and thus, we have shown that the PEPS with tensor ( 7.12 1, using a semi-regular 
representation, describes a state which is locally equivalent to the double model 
for the same group. 



8 Conclusions and outlook 

In this paper, we have presented a general framework for understanding Matrix 
Product States and Projected Entangled Pair states, and thus a large variety 
of interesting quantum states, in terms of symmetries. While this classification 
allowed us to reproduce the known results in one dimension, it also enabled us 
to go beyond one dimension and characterize the properties of two-dimensional 
systems with non-unique ground states. Using the symmetry-centered perspec- 
tive, we could show how these states arise as finitely degenerate ground states 
of local Hamiltonians, and characterize the different ground states using sym- 
metries at the closure. By looking only at these symmetries, we were able to 
explain the topological correction to the area law, the local indistinguishability 
of the ground state subspace, the anyonic nature of excitations of these systems, 
and why they are fixed points of renormalization group flows. 

We believe that the framework set forward in this article will allow to ex- 
plain a variety of things beyond what we have presented in this paper. For 
instance, the results on the parent Hamiltonian and ground state degeneracy 
can be directly extended to three and more spatial dimension (replacing the 
pair-conjugacy classes by triplet-conjugacy classes, etc.), as well as the state- 
ments made about the topological entropy, local indistinguishability, RG flows, 
and commuting Hamiltonians. By introducing excitations which are e.g. mem- 
branes of g's on the bonds, the framework of G-injective PEPS will also allow 
to explain the excitations of e.g. three-dimensional models. 

The results presented here likely can be generalized to models with other 
symmetries. For instance, on might think of tensors for which the ordered 
product of the group elements assigned to all indices is the identity for all non- 
zero entries. This symmetry shares the crucial property of being stable under 
concatenation, and it allows for the definition of injective and isometric tensors, 
and thus should allow for comparable statements about topological entropy, 
local indistinguishability, renormalization transformations, or non-abelian exci- 
tations. Note that Kitaev's code state and the double models can be expressed 
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in this form if one constructs the PEPS by starting from a ground state of Hb 
and subsequently projects onto the ground state subspace of Ha- More gen- 
erally, similar results might hold when considering more general symmetries, 
e.g. tensor categories used to define string-net models 
network based ansatzes 
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or different tensor- 
Note that from a different perspective, the 
role of so-called gauge-like symmetries for topological order has been studied 
in 
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The framework for understanding RG fixed points by mapping between iso- 
morphic representations suggests the extension to the case of renormalization 
flows, e.g. by considering the tensor product of individual bond symmetries. 
For instance, one could use the fact that the normalized character of the tensor 
product of any faithful representation converges to the character of the regular 
representation, up to global phases, to show that RG flows converge towards 
G-isometric PEPS. 
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